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The adiabatic, Holstein-Hubbard model describes electrons on a chain with step a interacting 
with themselves (with coupling U) and with a classical phonon field tfx (with coupling A). There is 
Peierls instability if the electronic ground state energy F[if) as a functional of tpx has a minimum 
which corresponds to a periodic function with period where pp is the Fermi momentum. We 
consider ^ irrational so that the CDW is incommensurate with the chain. We prove in a rigorous 
way in the spinless case, when A, U are small and y large, that a)when the electronic interaction 
is attractive U < Q there is no Peierls instability b)when the interaction is repulsive f/ > there 
is Peierls instability in the sense that our convergent expansion for F{(p), truncated at the second 
order, has a minimum which corresponds to an analytical and ^ periodic i^^- Such a minimum is 
found solving an infinite set of coupled self-consistent equations, one for each of the infinite Fourier 
modes of tp^- 



I. INTRODUCTION 

In 1955 Peierls, in Ref. suggested that in a one dimensional metal it is energetically favorable to develop a 
periodic distortion of the linear lattice with period ^ where is the Fermi momentum of the conduction electrons. 
The attempt of the conduction electrons to screen the periodic potential generated by the periodic lattice distortion 
creates a Charge Density Wave (CDW) in the conduction electron density. If a is the step of the undistorted lattice, 
then, depending whether 2^ is a rational number or not, the CDW (or the periodic lattice distortion) can be either 
commensurate or incommensurate with the non distorted lattice. While a commensurate CDW has preferred positions 
in the lattice, an incommensurate CDW has not and so it can slide without any change of energy; this was considered 
by Frolich in Ref. JTst a possible mechanism for superconductivity or a least for an enhancement of conductivity (see 
Ref. H). 

Indeed starting from the 70 's both commensurate or incommensurate CDW with wavevector 2pp have indeed been 
observed in a number of compounds (see for instance Refs. |3C| ] and fl^ ). A new wind of interest followed recently 
the discovery of high superconductors showing one dimensional incommensurate CDW, see Refs. Q and | |2^ . 

From a theoretical point of view, Peierls instability can be studied in the Holstein-Hubbard model, which is the 
simplest model involving both an electron-phonon and an electron-electron interaction. The standard theory of CDW 
is usually developed within the adiabatic approximation, in which the phonon field is treated as a classical field and 
the model becomes variational. We will consider the spinless case, so that the spinless adiabatic Holstein-Hubbard 
hamiltonian is 

H = Ho + Hp + XP + UV= (1) 
{txy - l^Sx,y) i'ti'y +\'Y'pI - Vxi^ti^x + U Y^'^^'^^ ~ \\\^t+i^x+i - 

x,i/eA xGA x^K kGA 

where A = 1, i — 1, tx^y — Sx,y — ^{5x.y+i + Sx,y-i), S''^^ fermionic creation or annihilation fields with periodic 
boundary conditions, is the chemical potential fi — 1 ~ cospp- In (|^) Hq is the fermionic kinetic energy. Hp is the 
phonon kinetic energy, XP is the electron-phonon interaction and UV is the electron-electron interaction; to describe 
the Coulomb repulsion one needs U > but it is not irrealistic to consider also f/ < (in this case UV is an effective 
interaction taking into account phonon- mediated processes). If C/ = the above model is called Holstein model, pp 
is the Fermi momentum of the non interacting U = model; the Fermi momentum in the U case is in general 
different, but we fix it to pp by adding a term i^N to the hamiltonian, where N — 'Ylxeh'^t'^x is the total particle 
number operator and u is a suitable counterterm. The proof of Peierls instability consists, within this model, in 

2 

the proof that the ground state energy F{lp) — X^^eA 2 ~^ ^oi'f) is minimized by ipx = (pl2ppx) where (p(t) is a 
27r-periodic function. The existence of a global minimum of the form (p{2ppx) was proved in the half filled band case 
Pf = tt/2 in the Holstein model and in the spinning Holstein-Hubbard model (see Refs.j20| and [0). Local minima 
of the form ip{2ppx), for any pp = tt^ with P, Q relatively prime, were found in Ref. for |A| < OilogQ^^). We 
are interested here in the incommensurate case in which ^ is an irrational number. At finite L, it is not possible 
to fix ^ directly to an irrational number, as in this way Cp(2ppx) cannot verify periodic boundary conditions. We 
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look however to a sequence of Li, rii such that liiiii^oo Li = oo and Unij^oo Pf.i = Pf and ^ is irrational, where 

Li 



PF,i = ■ An incommensurate phonon field has the form 



[(i.-l)/2] 

= lim y ^„e'2p,,.„.^ (2) 



n^-[L,/2] 



We require moreover that 



\2npFM\T^ > Co\n\-\ O^neZ, l"l < (3) 



where ||fc||xi = niin„g2 1^ ~ 27rn|. (|^) means that pf/t^ is a Diophantine number (see for instance Ref. p5|) and the 
proof of the existence oi pF,i verifying can be found in Ref. Such assumption is not really restrictive, as if 
T > 1 the complementary set of such points has measure 0. A Diophantine condition like (^) (in the infinite volume 
limit) appears in classical mechanics, for instance in the KAM theorem, see Refs. and |2^], and it is useful for 
handling with the so called small divisor problem; a similar problem appears also here. That a sort of extension of 
KAM-techniques to quantum systems is necessary to prove Peierls instability in the incommensurate case was pointed 
out for the first time by Aubry in Ref. |^ by analogy with the Frenkel-Kontorova models (see also Ref. ) . We have 
to specify the space of functions on which F{(p) is defined as a variational form. We say that F{ip) is a functional 
F : ^} —> R where fl is the set of functions (p^ of the form with zero average tpo = and ipn = ip-n = 'fn- 
Moreover, if k, Fq, Fi > are constants and a — X(pi then |cr| < Fq and, for \n\ > 1 

|A^„| <Fi|a|e-''l"l. (4) 

The condition (|) ensures also that the 2pp harmonic is present ( if (ySi = then (p^ is a constant). If (p is an 
extremal point of F{(p), it must satisfy the condition (po — Xp where p is the fermionic density. On the other hand, we 
can always include ipQ in the chemical potential p and then we can restrict our search of local minima of the ground 
state energy F{tp) to fields (p with zero mean. Note that F : U ^ R is indeed a function of the Fourier coefficients 
tpn] this means that, at finite L, can be considered not a functional but a L-dimensional function and only at the end 
we will take the L —^ oo limit. 

We say that there is Peierls instability if the variational form F : —^ R has a minimum ip^ €z ^ . We will show 
in §2 that if 1^9 j; e 17 is a local minimum then 

cpn = Xpn{^) n^O, n--[L/2],...,[(L-l)/2], (5) 

and Mnm = ^nm — Pm{'f) positivc definite, where Pn{^) = ■ Peierls instability can be proved by solving 

the infinite (as L — > 00) set of coupled equations (^). There are then two main steps to be performed; the first is to 
compute pn by an expansion, as there is no hope to compute it in a simple exact form, and the second is to solve the 
system ipn — Xpn{'p). In §2 an expansion for Pni^) is found, which is convergent for any ip G for A, U small enough. 
The proof of convergence is based on a sort of generalization of KAM theorem to quantum system; in §2 we review 
the main ideas referring to Ref. |Q for the mathematical proofs. The result is that we can write /5„ = J^kLo P-"^^ 
and, if pfa verifies and (p^ G fi, then for |A|, \U\ < e one has \p'k^\ < f{n,X,U)C^e^, where C is a constant and 
f{n, A, U) is a proper function. Of course the exact form of f{n, A, U) is important, and it will be specified in §2.E. 
The proof of the convergence (for e small enough) of the expansion for the density is based on Renormalization Group 
methods; it is important however to stress that, while the usual RG methods are only approximative, our results are 
mathematically exact; we refer for an introduction to rigorous RG for fermions to Ref. ||l6|] . 

The second step (see §3) consists in solving (j^). This is quite a difficult task, and we consider the simpler problem 
obtained by keeping only the first two terms of the expansion for p„. This means that we study 

= Ap(°)(^) + Ap(i)(^), n = -[i/2], [{L - l)/2]. (6) 

so neglecting O(e^) terms in the r.h.s. of (p|). The convergence of the expansion for /5„ makes this approximation 
quite reasonable. Note also that pl"^ and /5„ are rather complex function of A, U (the expansion is not a power series 
in A, U). Our main result is that if is small enough, in the L — )■ 00 limit 

a) in the attractive U < case there are no solutions (p £ of ^ (and of as well); 

b) in the repulsive case U > there is a (px £ ^ solving ^ such that 
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with 7] = PiU + 0{U'^) is a critical index and a, A positive constants. This means that, if the electron-electron 
interaction is larger than the electron-phonon interaction, there is a dramatic dependence on the repulsive or attractive 
nature of the electron-electron interaction. In the attractive case there is no Peierls instability as there are no functions 
analytical and ^ periodic minimizing the ground state energy. In the repulsive case there is instead Peierls instability, 
in the sense that the energy, keeping the first two non trivial terms of its convergent expansion, has a minimum which 
corresponds to an incommensurate CDW. It is very reasonable that higher order terms do not change this result (see 
the considerations in §3.B). 

Our results improve preceding works on the subject in which a)the interaction among electrons was neglected; b)the 
analysis was restricted to the first order 

^„ = Ap(o)M, (8) 

where, if U — O.Xpl^^ip) ~ —aX^ipi log \Xipi\ if n = ±1 and zero otherwise. There is a major difference between (^) 
and (|6|); while admits trivially a solution A(£i„ = a6n,±i and a is obtained by a BCS like equation, a solution of 
the form A(^„ = a-dn.±i does not solve (^) or (|^). While a solution of (|) if [/ = (when the model reduces to the 
Holstein model) is trivial to find, we are not able to find a solution of (^) £ fl ii U — 0, as the method we use to find 
a solution of (^ when U > and |j| is large fails in that case. We have then no evidence of Peierls instability in the 
Holstein model and it is unclear if this is merely a technical problem or if there is J7c > such that Peierls instability 
holds only for U > Uc in the Holstein-Hubbard model. 

The paper is organized in the following way. In §11. A we describe an expansion for the ground state energy of the 
Hubbard-Holstein model (|^) for ip £ il. In order to prove the convergence (for small A, U) of the expansion for small 
A, U, one has to solve a small divisor problem, and this is discussed in §11. B. In §11. C it is studied the RG flow and 
in §11. D we prove that the renormalization of the Fermi momentum if A-independent. From the ground state energy 
expansion, is it easy to derive an expansion for the density, and this is done in §11. E. In §111 we prove that (^ has no 
solution (fi £ fl in the attractive case, if |^ is small enough, while in the repulsive case we find a solution (ys e 17 of 
@ by a contraction method. Finally in §IV we discuss some open problems, in particular for the U = case. 

II. RENORMALIZATION GROUP ANALYSIS 
A. Grassman integrals 

It is well known that Eo{ip) can be written as a Grassman integral ( we use the same symbol ip for field and 
Grassman variables with a traditional abuse of notation) 

Eo{^) = - hm log / P(d^)e-^^-^^-''^, (9) 



/3^oo LP ' 



where 
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and X — [xq, x) and x + 1 = [xq, x + 1). P{d^p) is a Grassmanian integration defined on monomials by the anticom- 
mutative Wick rule with propagator 

^^'^'^y^^jL^ ~iko ~cosk + cos PF ' ^^^^ 

k 

where k = (fco, k). (|^) has a well defined L, /3 — > cx) limit only if the counterterm ly is chosen in a suitable way as a 
function of the parameters appearing in the Hamiltonian so that the Fermi momentum is just pp. In order to find 
the minima of F{ip) we have to differentiate with respect to (pm so one has in principle to take into account the 
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possible dependence of v from i^„, which is in general very complicated. However we will show in §2.D that it is 
possible to choose v as independent from A and so from (^„. This is due to the fact that the chemical potential can 
be moved inside the gap opened by (fx without affecting any physical property, and we can use this freedom to fix 
v as independent of ipn- It follows that a necessary condition for ip^ € fl to he a local minimum for F{ip) is that it 
verifies fx = Xpx where px — lim^^oo,r^o -zS^'^IxjT-^XjO) and S^'^{x,t;x,0) is the Schwinger function defined by, 

if (f>i are Grassman variables and writing J dx = /J^ ^xeA^ 

S'^'^i^-^y) = -^log / P{di^)e-''^^^-I''^'^i^^+^~^^\=o, (12) 



and V — UV + XP + vN. (||) is obtained by the Fourier transform of (px — \px defining px = X]n=-|-i/2] e^™P^^yO„. 

The above Grassman integrals can be evaluated by Renormalization Group methods; we refer to Ref. [p^ for 
an introduction to the formalism we are using and to Rcf. [|4| for the mathematical proofs of the convergence of 
the expansion we are describing. We start by evaluating the partition fmiction J P{d7p)e^^'^^\ It is convenient to 
decompose the Grassman integration P{dip) into a product of independent integrations. Let be \k\ = -\//cq + ||fc|||,. 
We write 

gik) - fi{k)g{k) + (1 - fi{k))g{k) = (13) 

where /i(fc) = 1 — xik — Pf, ^o) ~ x{k + Pf, ^o) and x(fc', k) is a C°° function with compact support such that it is 
1 for |fc'| < ^ and for \k'\ > qq, where 7 > 1 and 7, ao are chosen so that x{k i PF,ko) are non vanishing only in 
two non overlapping regions around ±p f ■ We write k ^ k' + ujp i? , w = ± 1 and 

h— — oo uj—±l h— — oo 

where fh{k') — xil^^k') — xil^'^^^k') has support 0(7'') around ujpF- The integration of the ultraviolet 

integration, gives 



e-v<o)(.0(^o)) ^ J ^^^^^(i)^^_ve^(i)+^(.o,)^ ^^^^ 
where ip-^ = Sfc^-oo (denoting simply by by J dk) 

00 00 „ 2n 
V(0)(V;(<")) ^Y.T. dfcl...rffc2„4f"^"^.41°^"^"^2n,™(^l,---,fen)^(^^,fc.+ (16) 

1^^^ — r\ ^ ^ — 1 



n— 1 711—0 ^ 



where pF — {pf, 0), Ci = ± and the kernels W^ „^{ki, A:„; z) are C°° bounded functions such that ,„ = 

and \WU < C"2'»'^'^(2."/2-i) ifz = Max{\X\, \U\, By an explicit computation it follows that 1^4% = U + 0{U^), 

= 0{Ua) for m 7^ and VF" r?j. = f + 0{aU) for m ^ 0. V*-"-* is called effective potential at scale 0; note that it 
contains non local interactions between an arbitrary number of fcrmions. 

The study of the infrared integration is much more involved. Following Wilsonian Renormalization group methods, 
we have to identify the relevant, irrelevant and marginal interactions; this is done by a power counting argument and 
it turns out, as it is standard in fermionic one dimensional systems, that the interactions quadratic in the fields are 
relevant, the quartic are marginal and the other processes are irrelevant. However there are in this model L many 
different terms bilinear or quartic in the fields, depending on the value of m in (|l^), and so it seems that there are 
L different non irrelevant interactions to be taken into account, which seems an hopeless task as we are interested 
in the L oo limit; this problem will be solved by a improved power counting in which the Diophantine condition 
plays a crucial role. Note finally that the quadratic interaction have a non trivial flow, so it is necessary to change 
the fermionic integration at each step; in other words the model has an anomalous behaviour, due to the fact that the 
model is close to a Luttinger liquid. 

The integration is performed iteratively, setting Zq = 1, ctq = cr, in the following way: once that the fields 
...,ip'^'^~^^^ have been integrated we have 

FzJdV.(^''))e-^"''(^'/''""). (17) 
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Then, putting Ch{k') ^ = X]j=-oo fji^') ^^'^ a{k') = (cos A:' — 1) cospp, dq = sinpp: 

fc' '^^il 

exp f - V / dk'C,,ik')Z^\( - iko - a{k') + cvvo smk')^f';^+ V^^^Z - ^/.(fc') V-^rTf^t ^f^^Z ] }■ 



Note that after \h\ steps the integration is different to the initial one; there is a wave function renormahzation Z/,, and 
a mass term ah- Moreover the effective potential at scale h has the form 

VW(^(<^)) = J2T. dk[...d%^ n ^£itS.,./(E ^^(^^ + ^^^^^) + 2mp^^)W^2t,™(fc; + ^iPF, ..; {^}). (19) 



n— m— 

In order to integrate '0'^^^ we write V^*) as rve^) + 7^V(''^ with 7^ = 1 - £. The £ operation is defined to extract the 
non irrelevant terms in V*-''-'; it is easy to check from a power counting argument that the terms in V^''-' involving six 
or more fields are irrelevant, thus £ = on such terms. Moreover, we will define £ = on the addenda in ( p^ not 
verifying the condition 

(TiUJipp + 2mpF = mod. 2tt, (20) 

which means that we are considering irrelevant the terms such that the sum of momenta measured from the Fermi 
surface is not vanishing (but it can be arbitrary small, due to the irrationality of ^). At the moment this definition 
of £ is completely arbitrary; it will be clear in the next section where we will show that the terms non verifying ( [20| ) 
are indeed irrelevant (here is where the Diophantine condition plays a role). 

In conclusion the definition of £ is the following: 

1) If 2n > 4 then 

CWt,raiki,-) = 

2) If 2n = 4 then 

jC.Wl„,iki, ...) = 6,n,oSj24^^^^^^^^Wl„,iuJiPF, ...^UJiPF) (21) 

3) If 2n — 2, ui = UJ2 then 

^{^2,ni{K +^lPF,k2 +UJ2Pf) = (Sm,0 (t^lpV , ti^2pV ) 

+ujiE{k' + ujipF)dkWl^{uJipF, uj2Pf) + k°dkoWl^(uJipF,uj2PF)], (22) 

where E{k' + upp) — vquj sin fc' + (1 — cos k') cospp and the symbol dk, dko means discrete derivatives. 

4) If 2n = 2, wi = -CJ2 then 

^W^.mik'l + i^lPF,k2 + 1^2Pf) = Sra.,^^W^^^{LUiPF,UJ2PF)- (23) 

The Kronecker deltas in the r.h.s. of (pl|), (p2[) ([2^) ensure that £ = if (^0|) is not verified. 
We find 

£V(")(^) = + Fl^^^ + + + a^Fi^^^^ + (24) 



where 



' -\-UjpF k' — UjpF , — 



I ^k'+ujpF,uj^k'- 

f/^'') = y f dk' Mk')i;i^"^+ , (25) 

UJ = ±1'' 
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where i = i',C,Oi and = 1, /c = —iko and /„ = E{k' + copp); moreover uq = U{v{0) — v{2pp)) + 0{U'^), 
So = 0{UX), oo, zo — 0{U^), no — v + 0{U). Note that in CV^'^^ there are terms renormahzing mass and the wave 
function renormahzation and it is convenient to include them in the fermionic integration writing 

where P2^_j (d?/'^^'''') is defined as (d?/; eq([l8|) with Zh-i and ah-i replacing Zh.Uh, with 

Zh-i{k') - Zh{l + C^\k')zh); Zn-i{k')an-i{k') - Zh{an{k') + Cy^\k')sH. (27) 
Moreover V^^' = CV^^^ + (1 - 'C)V('') and 

£V('') = 7""/.^^i-'^ + (a. - Zh)F(^"^ + U,P(^") (28) 
The r.h.s of ( p6|) can be written as 

Pz._,(dV(^''-^)) / Pz._,(dV'^''^)e-^"''(^'^'"") (29) 



where Pz^.^id^J^-''^^^) and Pz^_^{d%l;^^'^) are given by (|8|) with replaced by Zft_i(0) = Zh-i and Ch(fc') 

replaced with Ch-i{k') and f^^ik') respectively, if 

A(fc') = ^.-i[^^-^^] (30) 

and tp^-^'^ replaced with ip(^'^~^'> and ■0^''^ respectively. Note that fh{k') is a compact support function, with support 
of width 0(7'') and far 0(7'') from the "singularity" i.e. Lopp- The Grassmanian integration Pzh-iidip^'^'') ti^s 
propagator 

given by 

^ ( ^t'-ik'jx-v) 1 ^ -ikQ-a{k')-Vf)Sva.k' ah-i{k') \ , 

Zh-J A-i(fc') ^^V <^h-i{k') -iko-a{k')+vosmk' J ^ ^"^^^ 

where Ah{k') — [—iko — a(fc')]^ — (uosinfc')^ — [cr/i-i (fc')]^. It is convenient to write decompose the propagator as 

aitiS-y) = aliix-y) + C^^Hx-y), (32) 

where 

9^U^ -y)-jrpX^ -zfco+LosLfc' + a(/fc')^'^^'^ ^^^^ 



and, for any integer > 1 



Moreover 



i„w (f_^)i<i^i 



(34) 



(35) 



Finally we rescale the fields so that 
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(36) 



where 

CVC-) (V^) = ^"ly^.Fi^''^ + ShFi^'^^ + UhFlf''^ , (37) 

and by definition 

Vh = „ nh; 5h = {ah - Zh); Uh = ( „ fuh- (38) 
We perform the integration 

1 Pz._.(dV^''))e-^"''(V^'^'"") =e-^"'"^'(V^V.<^''-')^ (39) 

where V*-''"^-' has the same form as V^^-* and the procedure can be iterated, as the insertion of (|3^) in (^9|) gives an 
expression hke ( |l7| ) with h—\ replacing h. The above procedure is iterated untiU a scale h* defined as the minimum 

h such that 7'' > \oh\ is reached. Then we will integrate directly the field ^ — ^Y^=-oo'^^^^ without splitting 

the corresponding integration in scales (as was done for h > h*). This can be done as {x — y) verifies the bound 
eq( p4h with h* replacing h i.e. it verifies the bound valid for a single scale; the reason is thatfor momenta larger than 
0(7^) the theory it is essentially a massless theory and for momenta smaller is a massive theory with mass 0(7'' ). 
We will call running coupling constants Vh = {Uh, 6h, i^h) and renormalization constants Zh,crh', their behaviour as 
a function of h can be found by an iterative equation called beta function. Note that the irrelevant terms are not 
neglected, contrary to what is done in the usual RG methods, which are only approximative and not mathematically 
exact. The expansion generated by our RG is instead exact in a mathematical sense and nothing is neglected (see 
Refs. [H and ^ for details). 



B. The small divisor problem 



We have considered irrelevant the terms involving two or four fermions in the effective potential not verifying 
(|^). Looking at (|l|) we see that each addend contributing to the effective potential describes the interactions of 2n 
fermions whose momenta measured from the Fermi surface verify 

2n 

^ CTifc^ = ^ UiUJiPF + 2mpF. (40) 

2—1 i 

Then ( |20| ) says simply that the non irrelevant terms are only the ones in which the sum of the momenta measured 
from the Fermi surface is vanishing modulo 27r. This condition seems very natural in the commensurate case i.e. 
when pf = tt^; in that case, for n = 1, 2 if the r.h.s. of ( pO| ) is non vanishing modulo 27r, then it is greater than O(^) 
so for Q not too big at least one fermion involved has a momentum far enough from the Fermi surface. However 
things are not so simple in incommensurate case; in such case forn== 1,2 the r.h.s. of (|l|) can be very small (modulo 
27r) for very large m; in other words there are terms in the effective potential which are dimensionally relevant or 
marginal involving fermions with momenta arbitrarily close to the Fermi surface and not verifying (|20|); for instance 
V'fe+pp i'4'k-pF+2mpF -1 with 2mpF + 2kTT ~ for a suitable k. Then in the incommensurate case it is not clear if 
the terms not verifying (^ ) are really irrelevant (this problem is often not seen in literature, see for intance Ref. 
pl[) . This problem, with an interely different language, is well known in classical or celestial mechanics as the small 
divisor problem, for instance in the KAM or Lindstedt series for invariant tori of an hamiltonian system close to an 
integrable one. It is pos sible to write such classical series in terms of Feynman graphs so that the similarity becomes 
very clear, see Refs. and [|^; the crucial difference is that such graphs have no loops, contrary to what happens 
here. Another remarkable case in which small divisors appear is in the study of Schroedinger equation with a quasi 
periodic potential (very related to our problem in the U = case); if a Diophantine condition is assumed on the 
period there are quasi-Bloch states if A is small, see Ref. jl^], for suitable values of the quasi momentum, while the 
eigenstates are localized {Anderson localization) for large A, see Ref. p9[ |. 

The fact that the contributions to the effective potential not verifying ( |20| ) are irrelevant in a RG sense means that 
the perturbative series as function of \k,5k,Vk and and are convergent in a neighborhood of the origin. 
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To give a complete mathematical proof of the above statement is not straightforward, as one has to use determinant 
bounds for the fermionic truncated expectation; one cannot simply prove that each Feynman graph admits a finite 
bound as the number of Feynman graphs at order n is 0(n!^), so we refer to Ref. [Q. However the key idea why the 
terms not verifying (|20| ) are irrelevant can be understood from an analysis based on Feynman graphs. 

Each W2n^^ admits an expansion in terms of Feynman diagrams defined in the following way. A fc-th order diagram 
contributing to W2n''m '^^^ obtained from k graph-elements representing the addenda in (|3^) or in TZV^^^ dl^ ) by 
pairing the half lines (bearing indices h, uj, k'). The unpaired lines are called external lines, and to each paired line we 
associate a propagator 5^' ci'iK) (0); integrating the product of these factors over all the momenta fc^ of the paired 
lines we obtain the value of the graph, if the expression is multiplied by a suitable sign to take into account the Fermi 
statistic. A maximal connected subset of lines with scales > hy is called cluster with scale hy, and denoted by v. An 
inclusion relation can be established between the clusters, in such a way that the innermost clusters are the clusters 
with the higher scale, and so on; see the picture for an example of graphs with its clusters, pictured as boxes including 
the paired lines. The half-lines (contracted or not contracted) are emerging by the end-points, associated to Vh or to 
the kernels of TZV'^; if to an end-point is associated Vk , the minimal cluster containing it has scale k. 



To each Feynman graph it is associated by the above rule a value, and so to each subdiagram associated to each 
cluster; moreover, to the lines external to a cluster is associated a momentum k' smaller in modulus than one flowing 
in the lines internal to a cluster. We denote by v each cluster, hy its scale {i.e. all the lines internal to the cluster v 
have scale < hy and at least one has scale hy, and the external are larger), by Py the indices of the external lines, by 
|P^| their number and by Nypp the sum of the momenta k of the external lines. If w' is the minimal cluster enclosing 
the cluster v, the 7?. 7^ 1 operation produces one or more k' associated with the external lines (which are bounded 
by 'y^^' ) and one or more derivatives on the propagators internal to the cluster (giving one or more extra 7~'M. At 
the end (see for instance Ref. [^) one gets the following bound (essentially found by power counting using (Q) for 
bounding the propagators) for a graph with k vertices and 2n external lines with value G^2n k 

i V 

where z„ = 1 if |P„| = 4 and Ny ^ 0; Zy I if \Py\ = 2 and Ny = ^2, wi = -CJ2; Zt, = 2 if \Py\ = 2, iV„ = 0; Zy = 

in all the remaining cases. The factor Zy is due to the action on TZ on each cluster v. In order to sum over all the 

l-P I 

possible assignments of labels hy we need that — 2 + ^-^^ + Zy > and this is true for all the clusters with two or four 
external lines not verifying (^0|) (except when \Py \ = 2 and uji = —UJ2 in which case it is zero). We have to improve 
the bound in the cases in which —2+ + Zy < 0. The first case we can consider is \Py\ — 2 and uJi — —012, Ny — uj2- 

A similar cluster can be produced only if a)there is a non diagonal propagator g\^\^ internal to the cluster v; b)if 
there is one or more points associated to (pn internal to v. In both cases this means that there is in the boimds an 
extra factor ak"f~'^ for some scale k > hy, see (p5|), and 



Wk\ h-k Wh\ ^ Kl i(^h 

k 

and the factor ^fhO^^^) allows us to sum over hy. 



a.|7-'= = H7'-'^^^^^*'"'^ (42) 
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It remains to discuss the clusters with |P„| = 2, 4 not verifying (|2C|); here is where the Diophantine condition conies 
in. Given a cluster v, if OiOJipp + 2Nypp ^ mod. 27r then 

\N.\>C[}—^-\P.\]. (43) 

In fact, by the compact support properties of the propagators and the Diophantine condition 

|P„| |P„| 
aol''-' > \\J2''^''^\\T>\\^N,pF + J2<J,LJ,pF\\T>Coi2\P,\ + \N,\)-\ (44) 
1=1 1=1 



The meaning of ( p3| ) is quite clear; remembering that the momenta of the external lines is 0(7'*"'), ( |4^ ) says that 
the external fields can have momenta very close to the Fermi surface only if 7V„ is very large. The corresponding 
contribution is then very small, for the exponential decay properties of ipn- We can define a depth Dy defined in the 
following way: if v does not contain any cluster Dy — 1 otherwise Dy = 1 + maXynDyii where v" are the clusters 
contained in v. It is easy to see that 

-TT -k\m\ T r T r k\N^\ 

X{Wn^<e^\{e^Jle-^^. (45) 

i i V 

Using ( ^3| ) and the fact that Dy < —hyi + 2 we get 



where J|* is restricted to clusters with two or four external lines with Ny 0; choosing 7^2 ^ > 1 we can associated 
to each of this clusters a factor 



exp[-fc7^2''-''+3] < 7^''^'' < j^ih^'-i^^) (47) 

and this extra factor allows us to sum over the scale indices. 

In the above discussion of small divisor problem we have used the exponential decay of ipn but such condition could 
be probably relaxed to a power law decay i.e. |<^n| < for some integer N i.e. relaxing the analyticity condition 
for ifx to a differentiability one. This is quite reasonable, by analogy with KAM theorem which is valid not only 
for analytic but also for differentiable perturbations, see Ref. |^6|] , but the analysis would be probably much more 
involved. 



C. The flow of the running coupHng constants 

The equations for the running coupling constants are, for h > h* 

vh-i = ivh + Gl Uh-i = Uh + Gl 



CTh-i = Oh + G'J bh-\ = 5h + G^ 



Zh 



= 1 + (48) 



It is convenient to split g'1^\ with i — ii,a,i' into 

G'l{fih,'^h, Oh] A*o, I'o, o-q) = G]'^{^ih, Vh] Mo, t'o) + G'l'^{^ih, Vh, Mo, J^o, o-q) (49) 

where we have spitted gi'^i as in (|3^) and G\'^ contains no non-diagonal propagators and only the part g^'^^ ^ of the 

diagonal propagators; moreover, there are no vertices with m 7^ 0; in G^''' are all the remaining contributions. It is 
easy to check that for i = fi,i>, for maxk>h\vk\ < e. 



\Gr\<C[^re\ (50) 
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This follows from the bound (^) for and from the fact that fi are momentum conserving terms. For i = cr, by 

\Gl''^{ph,Vh,a-h] ■■■] ^J-o,vo,<ya)\ < C\Uhah\- (51) 



symmetry reasons, G]'^ = and 



We decompose, for i = 



G]'''{nh, fh; A^o, I'd) = G\'^{^.h\ mo) + G^'ifJ^h, i^h; a^o, t'o) (52) 



where the first term in the r.h.s. of (q2|) is obtained putting i>k ~ 0, k > h in the l.h.s. It is easy to see, from the fact 

)lus 

\Gl^\fih;-;f^o)<Geji. (53) 



that gi^ljix; y) can be divided in a even part plus a correction smaller than a factor 7 '4 , that for maxk>h\vk\ < £ 



On the other hand 



|G;,^''(A.,;...;/i,,)| <C£^7^ (54) 



as one can prove using the exact solution of the Luttinger model, see Refs. §, [§,and ll7|. Moreover we have that, 
for i ^ v, ji 

\G\'^{p.h,i'h] ^Jl.o,vo)\ <Gvh\Uh\'^ (55) 

Finally by a second order computation one obtains 

Gl'^ = ahUhiPi + G^'], Gl''' = Ul[p2 + GW (56) 

with /3i,/32 non vanishing positive constants and \G];^\ < G\Uh\, and \Gl''^\ < G\Uh\- 

By using the above properties we can control the flow of the running coupling constants. In fact, if \i>k\ < 



constants ci, C2, C3, C4, G such that, if A, u are small enough and h > h*: 



Gelji + '-^] for any k > h* (which will be proved in the following section), it follows that there exist positive 



\Uh-i -U\< GU""^- e-^'^i^^" < < e-^^i^^''; g-'^^^i^''* < Z^,_, < e-^-^^'". (57) 



D. Determination of the counterterm u 

We show that it is possible to fix v to a A-independent value; more exactly we show that it is possible to choose i' 
as in the A = case so that Vk is small for any k > h* . In the A = case, ak — and there are no contribution to the 
effective potential with m 0; calling Vk^j^k the analogous of Vk^P^k^ we can write 

^/. = 7-"+M^+ E 7'-'Gi''=(i>,A)], (58) 

k=h+l 

where Gl'''{i',fi) = Gl''{i'k, P-k, P-o)- We choose 

1^ = - 7'"'Gi-^-(^,/i) (59) 



fc— — 00 



then 



-7"' E 7'-'Gi-^-(^,A) (60) 



and \Dh\ < Ge^* , as by (53) 
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7^" ^ j''-'\Gl'\9,i^)\<C'sj-'^ 7M'<C£7*. (61) 

k— — oo k— — oo 

For the model with A ^ 0, for h > h* 

1 

= ^ 7'-'[G^^'(^^,M,a)+Gi-^:.,M)]] (62) 

k=h+l 

and, inserting given by (^9|), 

1 1 
- = 7-"+'{ E 7'-'G^''=(i^,M,^)+ E 7'-'[Gi'n^,M)-Gi''^-(:^,/i)]} (63) 

k=h+l k=h+l 

We prove that, for h > h* , 

\vh-VH\<eG{^)\ <e(^)2. (64) 

The proof is done by induction, assuming that it holds for scales > /i + 1 and proving by eq(|63|) that it holds for scale 
h. Looking at the first sum in ( |63|) and using (|50|),(^7|) 

7-"- E 7^-^|G^'^.,M)l<Gi7'V E 7'-^(^)^ 

k=h+\ k=h+l 

< i^rC2e' E I'-'i-r < C,eY-,?- (65) 

^ kM T 

Finally we can write Gl'''{v,^) — G],'''{i',jl) = J2k>k ^k,kj '^ith 

D-k^k = Gl'''{vk, ^J.k■, i^fc, Ms; Afc+i; ''o, Ao) - Gl'^{vk,Hk] %, Afc; /l^+i; i>o, Ao) (66) 
and, by the inductive hypothesis and ( ^l| ) 

E|i?S,.l <GrGe^E7^^'="'^(^)'^^2Ge^(^)'' (67) 

fe>/£ fc>fc ^ ^ 



SO that the last sum in eq(|63|) is bounded by Ce(^) with G = 4C, and (|64|) is proved; from that equation it follows 
that for h> h* , then < Ge, so that it is possible to have the flow of i^h is bounded choosing a A independent v. 



E. Bounds for the density 

An expansion similar to the one we have seen for the effective potential can be defined for the Schwinger function 
and hence for p„. One obtains 

1 C30 1 

dkgiUk) + E E /^i.'^'i (68) 

h—h* k—1 h—h* 

A.k (h) 

with Pi _i being sum over Feynman graphs similar to the one for the effective potential in which h is the lowest scale 
of the propagators and k the number of points. There are two kinds of contributions. The ones in which there are no 
points associated to Xipn, which means that there is at least a non diagonal propagator; such terms are bounded by 
G'^s'^\afi\Z^^ . The remaining ones (in which there is at least an irrelevant point) arc bounded by G^e^"f~ ^ the factor 
7~ coming from (^l|). 

In the same way for n 7^ 1 p„ = Yl'kLi^\=h* P^'^'^^ ■ expansion of /5„ with |n| > 1 are contributing 

only graphs with at least one irrelevant point Ai^„' ; in fact it is not possible, by conservation of momenta, to get 
a contribution to p„ with |n| > 1 from graphs containing only [/, (5, vertices (taking into account that there are 
non diagonal propagators, the difference of the external momenta of such graphs can be at most 0, ±2pi?). Hence 

Ip^^^^I <e-«/2|"l|a|C"=e'=7*. 

This concludes the construction of a well defined algorithm for computing the l.h.s. of (5) with any prefixed 
precision, if A, U are small enough; in the next section we try to find a solution ip ^ Q, of (5) or (6) by a contraction 
method. 
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III. SOLUTION OF THE SELF CONSISTENCE EQUATION 



A. Contraction mapping 

We have seen in §2.E that p„ = 6n,±iPn^ + X^fcLo where pn^ is the first addend in (|68| ) and pl^'' is given by 
graphs with one point associated to an irrelevat term \(pn'- We keep in the r.h.s. of instead of the full expansion, 
just the terms with /c 0, 1, forgetting the others for the moment. 

We find natural, in order to find a solution of (^, to proceed in the following way, calling <!>„ = Xipn for \n\ > 1: 

1) We consider a = Xtpi as a parameter and we study <i>„ = X'^pn\x,U,a,^), looking for a solution (f>(A, C/, ct) in 
the class il. 

2) If a non trivial solution <i>(A, U, a) is found, we insert it in (6) with n = I looking for a solution a function of A, U. 
According to the above strategy, we have to solve, for \n\ > 1 

= -A2c«(a)(^„ + ApW(A,[/,a,$), (69) 

where 

c(i) (a) = / dk{g[]l {%') -g^^l {k' + 2np» + [dil i^') E -^^k' + 2npp + (1 - co)pp) (70) 
+ E ^£%{k')~9'tl{k' + '2npp-{l-u,)pp)+j2 J^^^-'-^^'^ E -^/-'L{k' + 2npp) 

h—h* h — h* h'—h'* 

+ E ^ffi'l^^fc') E -^;&' + 2npp)+ E ^^9[!.%{k') E ^gtU{k' + {2n + 2u:)p;p)\] 

h—h* h' — h^ h—h* h' — h^ 

and 

1 



rn 



E y 7 E E '^2m-^iWi-c.2+c^^2«A(^™ / dA?'5_^'^^^,(fc05iV^'(^' + (2m + w^^ (71) 

"2. '"a 

where g'^^\k) — g"'"'{k) — g^lj{k') ii k — k' + ujpp- We shall rewrite ( |69| ) as 

= ^pW(A,C/,a,$). (72) 

1 + X^Cn 

and we will look for a solution $(A, ?7, cr) of such equation by applying an iterative method. 

For a fixed L, $ is a finite sequence of L — 3 elements, which can be thought as a vector in R^^^ . We consider the 
space T — C^(i?^^^) of C^-functions of (pm = 2,3, — 1. We shall define a norm in T for any A,ct different 
from zero 



,1,1 r "iii°siAii , ,9$„,,i 

|$||^ = sup e^^[|ari|$„(a)| + I— ^|] 

|n|>l I OCT J 



(73) 



We shall also define 

B = {$G.F:||$||^<1}; 
The solutions of (fz^ ) can be seen as fixed points of the operator Ta : T T , defined by the equation: 

(l + A2cL'^(a)) 



In the sum over m in ( |7lD m ^ n so that there is at least a non-diagonal propagator, so that by (57), ii uj'i ~ — cji 
(say). 



12 



E ^} jdk'£t-.Ak')9^XU^k' + {2m + u:[~u:^)pF)\< (75) 

h,h' — h* h—h* 

Considering ci^-* (^Oj), it is easy to see that the integrals in the first four terms of ( [70| ) are bounded by constants, as 
there is at least a non diagonal propagator or an ultraviolet one; for the fifth integral the bound is 

E I / E i;;^!^' + < c, ^ ^ < (™) 

h—h* h'—h^ h—h* 

A similar bound is found also for the last integral of ( |70| ) 

A2 



Then from the above bounds, if jri is smah enough one has 



1 + X^dn 



(1) 



< 2\^. (77) 



We find a solution of ( [72| ) as the fixed point $ of Ta in B. Note that B is invariant under the action of Ta($), as by 
(El) 

|T,,„(<I>)|<2CA2[ sup ^^]|A|^H<i|A|^H. (78) 
fe=±l,±2 |A|~ ^ 

so that for A small enough 

l|TA(<i>)|| <1 (79) 

Moreover, by a similar argument 

||T,($)-T,($')||^<i|l$-$'||^, (80) 

so that Ta is a contraction. It is also evident that Ta(0) G B. Hence, by the contraction mapping principle, there 
is a unique fixed point $ of in which can be obtained as the limit of the sequence defined through the 
recurrence equation ^C^+i) = (<&*•'''''), with any initial condition e B. If we choose = 0, we get 

oo oo 

< ^ _ < ^ — < ii<i>wii^, (81) 

i=l i=\ 

which immediately imply that < 1. Then if is large enough there is a solution <i> of (69) which is a function 

of A, f/, (j; in the following section we insert it in (6) with n = 1 looking for a solution cr which is function of A, U . 

B. Determination of \(p\ 

We have now to insert $(A, [/, a) found above in (6) with n = 1, in order to find cr in a self consistent way. We can 
rewrite (^ as 

(^1 = -A2c(^)(a)^i + Ap^') (A, [/, a, $), (82) 

~Ah) (p-s 



where 



cW(.)^ / dk'{ E ^^Mi^+.9«(fc')3l^](fc' + 2p.) (83) 

h=h' ^ 
u=±l h=h- ^ 



h=h* 
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and more explicitly 



r 



X,U)^J2 / ^fc' = (84) 







E 



Zh cr fcg + sin"^ fc' + (1 - cos k'y 



where ct/j, verify ( |57|) , and |ci| < C where C is a constant. 
Then by (|4|) and (|68|) we obtain 

^[(\^)-V2 _ l][a-' + Uh{X,U,a)]+ XahiX,U,a) = 1, (85) 

with |/i|, I/2I < C and 771 = fi\U + ?7i, 772 = /3i[/ + ?72, |fy2| < CC/^, and C, a, A are positive constants. ( p5| ) is a 
non-BCS or anomalous self consistence equation describing a superconductor whose normal state is a Luttinger liquid; 
if t7 = it reduces to a standard BCS equation. If the electron-electron interaction is larger than the electron-phonon 
interaction,?. e. if is small, thus the result is very sensitive to the sign of U. 



a) In the attractive case U < there is no solution |(t| < 1 if is small enough. In fact all the terms in the l.h.s. 
of (|8^) can be taken arbitrary small if U, X and are small enough so there is no |cr| < Fq solving (|85|). Note that 
this is also true considering (5), i.e. there is no solution G il of (5) if is small enough for any ip fl. 

b) In the repulsive case C/ > for small enough there is a solution of the form (Q), as L ^ 00. Note that 
this expression for the gap is very different with respect to the BCS like form of the U = case; in particular it is 

as large as >le~^'~"'°^"~'| >> Ae^ . A similar expression for the gap appears in the interacting Kondo problem, 
see Ref. |2^] and in superconductors whose ground state is a Luttinger liquid, see Ref. |^^. It is easy to check that 
corresponding to such solution the Hessian is positive definite. 

Finally I discuss shortly how to solve the full (||) without truncating the r.h.s. at the second order. At finite 
L we can use a contraction method to find a fixed point of (|^). However the fact that for the higher order terms 
we have only bounds, and we cannot make use their exact expressions like was done for has the effect that we 
can only get a weaker bound for the decay of the solution i.e. a power law decay |A(^„| < j^pr instead than an 
exponential one. This was also what we could get in Ref. in the commensurate case. In the incommensurate case 
this forbids to take the L —> 00 limit, just because we have proved in Ref. the convergence of /)„ in that limit only 
if Xifn have an exponential decay; on the contrary in the commensurate case convergence holds with a weak decay 
condition and the limit can be taken. This problem seems however merely technical, and it could be solved proving 
convergence under power law decay (see discussion at the end of §2.B ); this would allow to prove Peierls instability 
in the incommensurate case without truncating the density expansion. 



IV. CONCLUSIONS 



In the perturbative regime of small U there is Peierls instability then for U > Uc.i and for U < Uc.2 there is not. It 
is an interesting open question whether Ud — Uc.2 — Uc or not, and if they are equal weather Uc is greater, smaller 
or equal to zero. We are at the moment not able to answer to this question; we can only say that when J7 = or [/ 
small with respect to X^ the iterative method for solving (^ is not working. Let us consider in fact the U — case. 

r— . -[g+al 

In that case (^5|) with n = 1 is a BCS like equation and gives |cr| — Ae >■'-' , where \g{U, X)\ < C\X\. However we are 
not able to find a solution to the equation for $; the contraction method used in the large y case here fails as there 
are integers n such that 2npp is very close to 2pF modulo 2tt (of course analogous considerations can be done if is 

close to —2pp) i.e. \\2npF — 2uipf\\t can very small. For such n, Cn'^ ( [70| ) computed aX U — {i.e. Zh — l,ah — cr) 
can be written as 

c(i) = c« + b + + o{X), (86) 

cr 

where 6 is a constant. As 1 + X^ci — O(A^) we have that 
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instead O(A^) as it is in the preceding section; so it is not clear how a contraction method could be applied (an 

explicit computation shows that times the coefficient of ipn+2 tends to 1 when a„ ^ 0) . This fact could be not 

simply a technical problem linked to tlie method. The idea underling the analysis of the preceding section, which is 
essentially the Peierls idea, is that the process involving the exchange of 2pF is the dominant one, and the others are 
corrections. However, in the incommensurate case for very large n, 2pF or 2npF are almost the same, due to Umklapp 
scattering, so it is not clear physically why the first harmonic is the dominant one. Technically this means that there 
are large n for which p„ ~ pi (see ( |8^ ) in which b is negligible with respect to ci which is log-diverging) so that the 
self-consistence equation (pn — Xpn for such n or for n — 1 seems similar, and it is not clear why |(^„| << as it 
should be in order ip^ to be analytic. This problem is absent in the commensurate case; in that case, see Ref. it 

is possible to see that \cn \ < ClogQ so that for Q < e i^i then ^^^q^. < |A| and the contraction method will work, 
so that Peierls instability is proved. 

The mechanism why the above difficulty is avoided in the Holstein-Hubbard model with U > and large, so 
that Peierls instability is proved in the incommensurate case, is that the harmonic with n = 1 has a non trivial flow, 
and becomes larger, while ipn for |n| > 2 has no flow; this is in a sense a consequence of the diophantine condition. 

The difficulty in finding a solution to (||) when [7 = could mean that there is Peierls instability only for U > Uc 
with Uc > 0, but of course a deeper analysis is necessary to conclude. Note that in the d = oo limit of the Holstein- 
Hubbard model an incommensurate CDW is found, see Ref. (l9), only for U > Uc with Uc non vanishing positive. 
On the other hand the existence of Peierls instability when [/ = is usually supported by the analogy of the Holstein 
model with the Frenkel-Kontorova chain in which the extremality condition for the energy gives a dynamical system 
known as standard map. For this model the existence of an incommensurate phonon field minimizing the energy 
is an application of KAM theorem in the perturbative regime (and in the strong coupling regime of the so-called 
Aubry-Mather theorem, in which Lp^ is not smooth but it has infinite many discontinuities). 



V. REFERENCES 



[1] V.I. Arnold, Russ. Math. Surveys 18,9-36 (1963) 

[2] S. Aubry, P.Y. Le Daeron: Phys. D 8, 381-422 (1983) 

[3] S. Aubry, G. Abramovici, J. Raimbault: J. Stat. Pliys. 67, 675-780 (1992) 

[4] A. Bianconi et al; Stripes and related phenomena, edited by Bianconi and Saini, Kluwer, New York (2000) 
[5] G.Benfatto,G. Gentile, V.Mastropietro. J. Stat. Phys. 89, 655-708 (1997) 
[6] G.Benfatto,G. Gentile, V.Mastropietro. J. Stat. Phys. 92, 1071-1113 (1998) 
[7] F.Bonetto, V.Mastropietro: Commun. Math. Phys., 172, 57-93 (1995) 
[8] F.Bonetto, V.Mastropietro: Phys. Rev. B, 56, 3, 1-43 (1996) 

[9] G. Benfatto, G. Gallavotti, A. Procacci, B. Scoppola: Comm. Math. Phys. 160, 93-171 (1994) 
[10] A.H. Castro-Neto,F. Guinea: Phys. Rev. Lett. 80,18, 4040-4043 (1998) 
[11] E.Dinaburg, Y.G. Sinai; Funct. Anal, and its AppL, 9,279-289 (1975) 
[12] L.H. Eliasson; Math. Phys. El. Journ. 2 (1998) 
[13] F H. Frolich: Proc. Roy. Soc, A223,296, (1955) 
[14] G. Gallavotti; Comm. Math. Phys. 164, 145-156 (1994) 
[15] G. Gallavotti, Elementary mechanics, Springer (1983) 
[16] G.Gentile, V.Mastropietro:P/i7/sics Reports, 352, 273-437 (2001) 
[17] G. Gentile, B. Scoppola: Comm. Math. Phys. 154, 153-179 (1993). 

[18] R.H. Friend: in Solitons and condensed matter physics, Ed A. R. Bishop and T. Schneder, Springer, (1978) 

[19] J.K. Freericks, M. Jarrell cond mat/9502098 

[20] T. Kennedy, E.H. Lieb: Phys. Rev. Lett. 59, 1309-1312 (1987). 

[21] E.H. Lieb, B. Nachtergale; Phys. Rev. B,51,8,4777-4791 (1995) 

[22] D.Lee, J.Toner, Phys. Rev. Lett. 69 (23) 3378-3381 (1992) 

[23] P.A.Lee, T.M. Rice, P.W.Anderson; Solid State Comm. 14,703 (1974) 

[24] V.Mastropietro Commun. Math. Phys. 201 81-115 (1999) 

[25] V.Mastropietro Mod. Phys. Lett. 201 81-115 (1999) 

[26] J. Moser. Nachr.Akad. Wiss. Cottmgen Math-Phys 1-20 (1962)) 



15 



[27] R.E. Peierls: in Quantum Theory of solids , O.U.P., (1955) 

[28] A.Perali, C.CasteIlani,C. Di Castro,M.Grilli, Phys. Rev. B 54,16216 (1996) 

[29] L.Pastur, A. Figotin: Spectra of Random and quasi-periodic operators, Springer, Berlin (1991) 

[30] G.A.Toombs: Phys. Rep. 40C,181 (1 978) 

[31] J.Vidal, D.Moihanna,T. Giamarchi: |cond-mat/990508C 



16 



